Experimental Limit on the Cosmic Diffuse Ultra-high Energy Neutrino Flux 
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We report results from 120 hours of livetime with the Goldstone Lunar Ultra-high energy neutrino Exper- 
iment (GLUE). The experiment searches for < 10 ns microwave pulses from the lunar regolith, appearing in 
coincidence at two large radio telescopes separated by 22 km and linked by optical fiber. Such pulses would 
arise from subsurface electromagnetic cascades induced by interactions of > 100 EeV neutrinos in the lunar re- 
golith. No candidates are yet seen, and the implied limits constrain several current models for ultra-high energy 
neutrino fluxes. 



In 1962, G. Askaryan predicted that electromagnetic cas- 
cades in dense media should produce strong coherent pulses 
of microwave Cherenkov radiation 1 1 ] . Recent confirmation 
of this hypothesis at accelerators I01 strengthens the motiva- 
tion to search for such emission from cascades induced by 
predicted high energy neutrino fluxes, closely related to the 
measured fluence of ~ 10 20 eV cosmic rays in many models. 

Two such models, the Z-burst model 1 3], and a generic class 
known as Topological Defect (TD) models @], predict ultra- 
high energy (UHE) neutrinos with either monoenergetic or 
very hard energy spectra. In the Z-burst model, UHE neu- 
trinos annihilate with relic cosmic background neutrinos via 
the vv — * Zo channel. The Zq then decays rapidly in a burst 
of hadronic secondaries which create the observed ~ 10 20 eV 
cosmic rays. The need to match the observed UHE cosmic ray 
fluxes and satisfy the current constraints on neutrino masses 
(which modify the annihilation resonance energy) then lead 
to a requirement on minimal neutrino fluxes at the resonance 
energy near 10 22 ~ 23 eV. The Z-burst model thus formally re- 
quires only neutrinos at a single energy, with no specification 
for how such a flux might be produced. 

The Z-burst model is also significant in that it is a variation 
on an earlier idea ||5[ in which the vv annihiliation process 
could be used as a probe of the cosmic background neutri- 
nos, one of the few viable ways ever proposed for detection 
of these relic cosmological neutrinos-it requires only a suffi- 
cient flux of UHE neutrinos and a detector with the sensitivity 
to measure them. Constraints on these UHE v fluxes thus can 
rule out this potential detection channel for the relics, in addi- 
tion to excluding their role in UHE cosmic ray production. 

TD models, in contrast, postulate a very massive relic par- 
ticle from the early universe which is decaying in the current 
epoch and producing secondaries observed as UHE cosmic 
rays. The required masses approach the Grand-Unified The- 
ory (GUT) scale at ~ 10 24 eV, and the decay products thus 
have a very hard spectrum extending up to the rest mass en- 
ergy of the particles. Because of these very hard spectra, 
detectors optimized for lower-energy neutrinos, even up to 
PeV energies, do not yet constrain these models, and new ap- 
proaches, such as the experiment we report on here, are re- 
quired. 

Neutrinos with energies above 100 EeV (1 EeV = 10 18 eV) 
can produce cascades in the upper 10 m of the lunar regolith 
resulting in pulses that are detectable at Earth by large radio 



telescopes, jfj @] One prior experiment has been reported, 
using the Parkes 64 m telescope |8] with 10 hours of live- 
time. In the decimeter band, the signal should appear as 
highly linearly-polarized, band-limited electromagnetic im- 
pulses fSJ]. However, since there are many anthropogenic 
sources of impulsive radio-frequency interference (RFI), the 
primary problem in detecting neutrinos is to reject such in- 
terference. Since 1999 we have conducted a series of exper- 
iments in search of such pulses, using the JPL/NASA Deep 
Space Network antennas at Goldstone, CA fioil . We have es- 
sentially eliminated RFI background by employing two anten- 
nas in coincidence. 

Although the total livetime accumulated in our experiment 
is a relatively small fraction of what is possible with a dedi- 
cated system, the volume of material to which we are sensitive 
is enormous, exceeding 100,000 km 3 at the highest energies. 
The resulting sensitivity is enough to begin constraining some 
models for diffuse neutrino fluxes at energies above 10 20 eV. 
We report here on results from 120 hours of livetime. 

The lunar regolith is an aggregate layer of fine particles and 
small rocks, consisting mostly of silicates and related miner- 
als, with meteoritic iron and titanium compounds at an av- 
erage level of several per cent, and traces of meteoritic car- 
bon. Its depths are 10 to 20 m in the maria and valleys, but 
may be hundreds of meters in portions of the highlands 1 1 1 ] . 
It has a mean dielectric constant of e = n 2 ~ 3 and a den- 
sity of p ~ 1.7 g cm~ 3 , both increasing slowly with depth. 
Measured values for the loss tangent vary widely depending 
on iron and titanium content, but a mean value at high fre- 
quencies is tan 8 ~ 0.003, implying a field attenuation length 
L a ~ 9 mat 2 GHz d. 

Fig- HI illustrates the signal emission geometry. At 100 
EeV the interaction length of a neutrino is about 60 km 
Upon interaction, a cascade ~ 10 m long forms, and Comp- 
ton scattering, positron annihilation, and other processes lead 
to a ~ 20% negative charge excess. This cascade radiates a 
cone of coherent Cherenkov emission at an angle from the 
shower axis of 9c = cos -1 (l/p^/e) ~ 54°, with an angular 
spread of AG ~ 1° at 2 GHz. At an ideal smooth regolith sur- 
face, the refraction obeys Snell's law, and the exit angle 0' 

is magnified by the gradient A0'/A8 = «cos8/ y 1 — « 2 sin 2 9 
which equals n for normal incidence (9 = 0), but becomes 
much larger as 9 approaches the total-internal reflectance an- 
gle Q T i R = ji/2— 9c- This magnification improves the accep- 
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FIG. 1: Geometry of lunar neutrino cascade event detection. 



tance for an earth-based detector, at the expense of an increase 
in energy threshold. Our ray tracing shows that similar effects 
obtain statistically for a more realistic surface as well. 

For our search we used the shaped-Cassegrainian 70 m an- 
tenna DSS14, and the 34 m beam-waveguide antenna DSS13, 
separated by 22 km. The S-band (2.2 GHz) right-circular- 
polarization (RCP) signal from DSS13 is filtered to 150 MHz 
bandwidth and down-converted to an intermediate frequency 
(IF) near 300 MHz. The band is subdivided into high and low 
frequency halves with no overlap. The DSS14 dual polariza- 
tion S-band signals are down-converted to the same 300 MHz 
IF, and a combination of bandwidths from 40-150 MHz are 
used for sub-band triggering on impulsive signals. At DSS14, 
an L-band (1.8 GHz) feed which is off -pointed by ~ 0.5° pro- 
duces a 40 MHz bandwith monitor of terrestrial interference 
signals. 

Fig-Elshows the layout of the trigger. The signals from the 
two antennas are converted to unipolar pulses using tunnel- 
diode detectors with a ~ 10 ns integration time. A comparator 
then test for pulses above threshold, and a local coincidence 
within 50 ns is formed among the channels at each antenna. 
The DSS14 coincidence between both circular polarizations 
ensures that the signals are highly linearly polarized, and the 
split-channel coincidences ensure that the signal is broadband. 

A global trigger is formed between the local coincidences 
of the two antennas within a 150 /is window, which en- 
compasses the possible geometric delay range for the Moon 
throughout the year. Although use of a smaller window is 
possible, a tighter coincidence is applied offline and the out- 
of-time events provide a large background sample. Upon the 
global coincidence, a 250 /js record, sampled at 1 Gsamples/s, 
is stored. The average trigger rate, due primarily to random 
coincidences of thermal noise fluctuations, is 3 x 10~ 3 Hz. 
Terrestrial interference triggers are on average a few percent 
of the total, but we have averaged > 95% livetime during the 
runs to date. 

The precise geometry of the experiment is a crucial discrim- 
inator for events from the Moon. The relative delay between 
the two antennas is x = c _1 |B|cos9 where is the apparent 
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FIG. 2: GLUE trigger and data recording system. 



angle of the Moon with respect to the baseline vector, B. For 
our 22 km baseline, we have a maximum delay difference of 
Xmax = ±73 /js. Detectable events can occur anywhere on the 
Moon's surface within the antenna beam. This produces a 
possible spread of 630 ns in the differential delay of the re- 
ceived pulses at the two antennas. 

The 2.2 GHz antenna beamwidths between the measured 
first Airy nulls are 0.27° for the 70 m, and 0.56° for the 34 m. 
We took data in three configurations: pointing at the limb, the 
center, and halfway between. The measured source tempera- 
tures varied from 70K at the limb, to 160K at the Moon center, 
with system temperatures of 30-40K. 

Timing and amplitude calibration are accomplished in sev- 
eral steps. We internally calibrate the back-end trigger system 
using a synthesized IF pulse signal, giving precision of or- 
der 1 ns. We use a pulse transmitter (single-cycle at 2.2GHz) 
aimed at the antennas to calibrate the cross-channel delays of 
each antenna to a precision of 1 ns. The cross-polarization 
timing at DSS 14 is also checked since the thermal radiation 
from the limb of the Moon is significantly linearly polarized 
(from differential Fresnel effects I14lll5l0 . introducing an eas- 
ily detectable LCP-to-RCP correlation. 

Dual antenna timing calibration is accomplished by cross- 
correlating a 250 /js thermal noise sample of a bright quasar, 
typically 3C273, recorded from both antennas at the same 
time and in the same polarization, using the identical data 
acquisition system used for the pulse detection. This proce- 
dure establishes the global ^136 /js delay between the two 
antennas to better than 10 ns. Amplitude calibration is ac- 
complished by referencing to a known thermal noise source. 
The system temperature during a run fixes the value of the 
noise level and therefore the energy threshold. We also check 
the system linearity using pulse generators to ensure that the 
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FIG. 3: An example of a GLUE triggered event. 

entire system has the dynamic range to see large pulses. 

Figure shows a typical event which triggered the system. 
The top two panes contain the DSS14 LCP and RCP signals, 
and a narrow pulse is present in both polarizations, indicat- 
ing a broadband spectral content, and a high degree of linear 
polarization. The pulse power is normalized to the local mean 
power over a 250 /js window. In the third and fourth panes, the 
two channels from DSS13 are shown. In the fifth pane from 
the top the L-band offset feed signal from DSS 14 is shown, 
and no RFI is present. The measured delay relative to the ex- 
pected Moon time is — 1 . 1 /js in this case, slightly larger than 
allowed. Systematic timing offsets from channel to channel 
are well under 10 ns. 

Two largely independent analyses look for pulses corre- 
sponding to an electric field 6o above thermal noise in all 
channels. Both analyses remove terrestrial RFI with either 
a visual or algorithmic method. Each enforces ~ 20 ns local 
coincidence timing at each antenna. The precise values of the 
cuts were determined and fixed before looking at more than 
half of the data. We have seen no candidates in 120 hours of 
livetime. The signal efficiency of the RFI and timing cuts is 
estimated to be > 98%. To estimate background levels, we 
also search of order 100 different delay values that are incon- 
sistent with the Moon's sky position, and we have also found 
no candidates in this search. Hence we observe no events with 
a background of < 0.01 events. 

The field strength (V/m/MHz) from a cascade of total en- 
ergy Wt in regolith material can be expressed |9]: 

1.45 xl0~ 7 W T v 1 
~ R TfeV Vo" (l + Cv/vo) 1 - 44 ) ° 

where R is the distance to the source in meters, v is the radio 
frequency, and the decoherence frequency is Vq — 2500 MHz 



for the regolith (Vo scales mainly by radiation length). For 
typical parameters in our experiment, a 10 20 eV cascade will 
result in a peak field strength at Earth of E ~ 1.3 /jV m 
for a 70 MHz bandwidth. We have verified Equation \l\ to 
within a factor of 2 through accelerator tests |2[ using sil- 
ica sand targets and y-ray-induced cascades with composite 
W T ~ 10 19 eV. 

Based on the effective antenna aperture r|A and the fact that 
the background events are due to fluctuations in the black- 
body power of the Moon, we estimate that the minimum de- 
tectable field strength for a linearly polarized pulse is 



I 2kT sys Zo 



(2) 



where Zq = 377£2, and N a the number of standard deviations 
required per channel relative to thermal fluctuations. For the 
lunar observations on the limb, which make up about 85% 
of the data reported here, T sy!i ~ 110 K (including the source 
contribution), v = 2.2 GHz, and the average Av ~ 70 MHz. 
For the 70 m antenna, with r| ~ 0.8, the minimum detectable 
field strength at A^j = 6 is E m i n ~ 0.8/jV m _1 . The estimated 
cascade threshold energy for these parameters is W,/„ = 6 x 
10 19 eV. Since the mean inelasticity is (y) — 0.2, the threshold 
neutrino energy is ~ 3 x 10 20 eV. 

Effective volume and acceptance vs. neutrino energy has 
been estimated via two independent Monte Carlo simulations, 
including the current estimates of both charged and neutral 
current cross sections 11311 . and the y distribution. The neu- 
trino species are assumed to be fully mixed upon arrival. All 
neutrino flavors were included, and LPM effects in the shower 
formation were estimated |9]. At each neutrino energy, a dis- 
tribution of cascade angles and depths with respect to the local 
surface was obtained, and a refraction propagation of the pre- 
dicted Cherenkov angular distribution was made through the 
regolith surface, including absorption, reflection, and rough- 
ness effects, and thermal noise fluctuations in the detector. 

Via ray-tracing in our simulations (see Ref. 11611 ). we find 
that, although the specific flux density of the events are low- 
ered by a factor of order 2 from refraction and scattering, the 
effective volume and acceptance solid angle are increased by 
as much as an order of magnitude. The neutrino acceptance 
solid angle is thus about a factor of 20 larger than the appar- 
ent solid angle of the Moon itself. UHE cosmic rays could 
produce a background for lunar neutrino detection; however, 
there are several processes that suppress cosmic ray radio 
emission with respect to that of neutrinos, including total in- 
ternal reflection, and formation zone effects II 611 . The level of 
the potential UHE cosmic ray background is unknown, but as 
we have detected no events of any kind, our limits obtain in 
any case. 

Figure 0] plots the predicted fluxes of ultra-high energy 
(UHE) neutrinos from AGN production lll7l Il8il . a maximal 
flux from UHE cosmic -ray interactions 1 19], two topological 
defect models ||4j , and the Z-burst scenario 1 3 ] . Two other cur- 
rent limits in this energy regime are plotted, from the RICE 
experiment |H, and the FORTE satellite JH. Our 90% 
confidence level, differential model-independent limit for 120 
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FIG. 4: Model neutrino fluxes & limits from the present work and 
other current experiments. The inset expands the region near the 
Zburst predictions. 



hours of livetime is shown plotted with large squares 1 23], 
based on the observation of no events above an equivalent 6o 
level amplitude (referenced to the 70 m antenna). 

All of the differential limit curves in Fig. |4] (except for 
the Ey^ power law AMANDA limit) correspond to the in- 
verse of the energy-dependent exposure = (neutrino aperture 
x time) for each detector noted, scaled by the Pois- 
son upper-limit factor (2.3 for 90% CL). Our plotted model- 
independent limit curves are very conservative, and become 
more restrictive when integrals over specific broad-spectrum 
models are considered. When converted to integral form, our 
limits constrain power-law neutrino models that parallel our 
curves for ~ 1 decade, even a factor of 3-5 below our plotted 
differential limits. 

For quasi-monoenergetic models such as the Z-burst sce- 
narios, the limits are exact, and the combined GLUE and 
FORTE results constrain much of the available Z-burst param- 
eter space. Models such as the highest TD curve shown 1 4] are 
largely excluded at higher energies by the GLUE and FORTE 
results, since the integral (flux x aperture) would lead to ~ 3 
and ~ 10 events in each detector, respectively, in conflict with 
the null results. 
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